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ABSTRACT 


Monthly  mean  vdlu<=s  of  temperature  output  for  a  ten-year 
period  (1969-78)  from  a  three-dimensional  primitive  equation 
ocean  model  are  used  to  compute  isotherm  displacements  at 
eight  selected  grid  points  in  the  North  Pacific  Ocean.  The 
ocean  model  input  parameters  are  climatological  radiation, 
clouds,  surface  air  tempeiature,  surface  humidity  and  six- 
hourly  synoptic  winds  taken  from  Fleet  Numerical 
Oceanography  Center  analyses.  The  firsr  five  baroclinic 
modes  are  computed  from  the  mean  temperature  profiles  at  the 
eight  selected  model  grid  points  and  these  modes  are  -hen 
fit    in  the    time    domain   to   the   isctherm   displacements. 

The  fit  shews  evidence  of  baroclinic  free  Rossby  waves 
in  the  time  domain  ar  the  selected  grid  points  south  of  40N 
and  indicates  the  predominance  of  the  first  baroclinic  mode. 
At  the  selected  grid  points  north  of  UON,  the  results  are 
inconclusive.  This  is  due  to  the  longer  periods  required  at 
the  northern  latitudes,  the  neglect  of  salinity,  and  masking 
by  a  strong  seasonal  signal  caused  by  surface  heating  and 
cooling   at    higher  latitudes. 
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I.        INTRODUCTION 

Oceanic  general  circulation  models  have  been  under 
development  for  the  past  decade  and  a  half.  Evaluation  of 
these  newly  developed  models  was  initially  hindered  by  a 
lack  of  adequate  observational  data  from  ocean  regions. 
However,  ever  the  past  eight  years  high  quality  oceanic  data 
and  data  analysis  studies  from  portions  of  the  ocean  have 
been  made  available  to  the  large-scale  ocean  modeller.  Of 
particular  imccrtance  and  interest  in  this  study  has  been 
the  analyses  cf  temperature  data  resulting  from  the  "North 
Pacific    Experiirent"     (NORPAX)     begun    during    the    last    decade. 

Eased  upon  an  analysis  of  TPANSPAC  data  (a  sub-program 
of  NCEPAX  confined  to  a  region  3  0-50N  and  170E  -  150W),  Kang 
and  Magaard  (1980)  determined  that  annual  temperature  fluc- 
tuations in  the  area  30-aON,  16 OE  -  130W  consisted  of  first 
mode  Rossby  waves  with  random  phase  traveling  northwest  with 
a  wavelength  cf  300  km.  For  other  latitudes  of  the  North 
Pacific  Ocean,  Price  and  Maqaard  (1980)  determined  that 
first  mode  baioclinic  Rossby  waves  with  dominant  period  of 
seven   to      ten   years    existed  in      the   region   40-50N.  It    was 

also  determined  that  Rossty  waves  of  annual  period  did  not 
play    a   role   in   the  latitude  belt    a0-50N  in   the    North    Pacific 
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Ocean.  Magaard  (1983),  in  a  paper  discussing  baroclinic 
Rossfcy  wav€  energetics,  summarizes  these  earlier  studies  by 
stating  that  first  mode  baroclinic  Rossby  waves  dominate  the 
internal  temperature  fluctuations  in  the  region  20-30N  in 
the  North  Pacific  Ocean  fcr  a  whole  range  of  periods  between 
five  months  and  ten  years.  Between  30-UON  only  annual  or 
near-annual  first  mode  baroclinic  Rossby  waves  are  evident. 
It  is  the  results  of  these  studies  that  have  ractivated  this 
investigation,  which  is  essentually  an  evaluation  of  whether 
a  general  circulation  ocean  model,  driven  by  observed  synop- 
tic winds,  can  simulate  the  above  characteristics  of  the 
North  Pacific  Ocean  Rossby  wave  field.  Following  a  proce- 
dure outlined  by  Emery  and  Magaard  (1976) ,  this  study  exam- 
ines the  possible  existence  of  baroclinic  Rossby  waves 
inferred  frcm  temperature  fields  produced  by  the  ocean 
model. 
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II.        MCCEL   DESCRIPTION    AND    ANALYSIS     PROCEDURE 

A.        MODEL    DESCRIPTION 

The  medal  used  in  this  study  is  a  three-dimensional 
primitive  equaTiicn  model  developed  at  the  Naval  Postgraduate 
School  and  described  by  Haney  (1980) .  A  flat  bot-om  and 
rigid  lid  are  frescribed.  The  model  is  driven  by  seasonally 
varying  (monthly  climatology)  solar  radiation,  clouds,  sur- 
face air  temperature,  and  surface  relative  humidity.  The 
surface  winds  utilized  are  taken  from  six-hourly  synoptic 
analyses,  1969-1978,         cade        at        the      Fleet        Numerical 

Oceanography   Center,    Monterey,    California.      Salinity    is    neg- 
lected.     The    model  domain    shown    in    Figure    1    is   a   rectangular 
basin   extending      from   the      Equator   to      65N   and      from    145E 
125W.  The    horizontal      grid   spacing      is    approximately      two 

degrees  in  latitude  and  three  degrees  in  longitude.  The 
model  is  4000  m  deep  with  20  irregularly  spaced  levels,  14 
of  which  are  above  1000  m.  The  boundary  conditions  pre- 
scribed at  the  bottom  and  sides  are  those  of  an  enclosed 
thermally    insulated    basin. 
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Horizontal  doiain  of  temperature  field  used  for  analysis  is 
shown  by  enclosed  region.  Triangles  mark  grid  points 
selected      for      Rcssby   wave  analysis. 

Figure   1.       Horizontal   Domain   of    Ocean    Model 


The  mcdel  also  includes  a  dynamic  adjustment  mechanism 
as  described  by  Adamec,  et.  al.  (1981).  Essentially,  the 
mechanism  assures  that  if  the  local  gradient  Richardson  num- 
ber. Hi,  falls  below  a  critical  value  (0.25)  rhen  the  temp- 
erature and  the  velccity  components  at  uhe  two  adjacent 
levels  are      mixed  vertically   so      that    heat    and      momen-um    are 
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conserved,        the    mixing      ratios    for      heat      and    momentum      are 
equal,   and   after   mixing    Ri   equals  the    critical   value. 

The  model  was  run  for  a  30-year  period,  with  the  first 
two  decades  considered  as  a  "spin-up"  period  and  the  third 
decade,  1969-78,  considered  as  a  simulation  of  actual  condi- 
tions. The  purpose  of  this  thesis  is  to  analyse  the  temper- 
ature variability  simulated  by  the  model  during  this  last 
decade.  The  ccal  is  to  see  if  the  variability  is  consistent 
with   the    dispersion    relation    for    baroclinic   Rossby   waves. 

B.        ANALYSTS    EBOCEDaRE 

'' •      Governing  Equaticns  for    Baroclinic    Ross b^   Waves 

The  governing  equations  for  baroclinic  Rossby  waves 
are  derived  from  the  frictionless  momentum  equations,  the 
hydrostatic  balance,  and  the  continuity  equation  for  an 
incompressible   fluid    in    the  following    form: 


u^    +    (V-V)u   -    ^v    +    p    Px    ""    °    ' 


v^    +    (V-V)v    +    fu    +    -  Py    =    0     , 


i  p      =  0    ,  (la) 

i  p      =  0    ,  -  (lb) 

V»V    =  0     ,  (Ic) 

P^    "^    PS    =  0     » 

p^    +    (V-V)p    =  0     , 


(Id) 
(le) 
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where  V  is  th€  velocity  vector*  u,  v  and  w  ar^  the  velocity 
components  in  the  eastward  (x) ,  northward  (y) ,  and  upward 
(z )  direction,  respectively,  p  is  the  density,  p  is  the 
pressure,  g  is  the  acceleration  due  to  gravity,  f  is  the 
Coriclis  parameter,  t  is  time  and  V=  i^kfSi'&z^'  '^^-  sub- 
scripts denote  partial  differentiation  with  respect  to  the 
subscripted  variable. 

The  vcrticity  equation  is  formed  using  equations 
(1a)  and  (lb).  Subsequent  application  of  the  quasi- 
geostrophic  approximation  and  linearizarion  leads  to  the 
follcwing  form  of  the  vorticity  equation: 


vp^^  +  p.'--) 


XX 


yy  t 


ep'  -  p  f  w'   -  0 


(2) 


where  p'  and  w»  are  the  perturbation  pressure  and  vertical 
velocity  ccraccnent,  0  =  df/dy,  p,  is  a  constant  reference 
density,  and  the  other  variables  are  as  defined  before. 
Additional  assumptions  used  in  the  above  linearization  pro- 
cess are  that  the  lateral  shear  of  the  basic  flow  is  neg- 
lected and  the  basic  staxe  horizontal  velocity  is  at  rest 
(i.e.,  u  =  V  =  0) .  From  this  assumption,  equation  (1c)  and 
the    boundary   condition   that   w   =    0   at    z   =   0    (sea    surface)       it 
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follows  that  w  =  0.  Equation  (2)  contains  two  unknown  vari- 
ables, p'  and  V*,  therefore,  an  additional  equation  is 
required. 

Substitutinq  from  (Id)  into  {1e)  after  multiplyinq 
by  "g/Po  si^d  subsequent  linearization  yields  the  second 
equation. 


p'      +  p    N^w'    =0  (3a) 

^zt      ^o 


or 


"'  =  -  ^^^  Pit  '  "" 

N    p 
o 


where  N"^  =  -(g/po)Pz  is  the  Brun t-Vaisala  frequency  squared, 
assumed  to  be  a  function  of  z  only.  Eliminating  w'  in  (2) 
using    (3b)    yields, 

N 

Using  the    separation    of    variables   technique   the    perturbation 
pressure    is    written. 


p'(x,y,z,t)    =    P^(z)    p^(x,y,t)     ,  (5) 
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where   n    is    the    mcde  number.      Substituting     (5)    into   eqaa-icns 
(4)     and    (3fc)    gives, 

1      d  r  1  n,     _  ^n    t "^n    x    _    .     I  /g) 

n           N  f    (Pj^)^ 


and 


1      dPn         .  Po    "' 


>,2    dz  f '^    \ 


(7) 


2  =2         -,2  2 

where      V  =    (^2*^2)  '      ^^<^    '^n  --    ^      ssparation   constant    with 
dimensions    cf      inverse   velocity.  Differentiating  equation 

(6)     with    respect    to    z   yields, 

,2 


d 


[1      dPn^     .     X     2    ^"=     0     .  (8) 


,    2    ^.Z    dz'  n      dz 

dz        I1 


Substituting'  frcoj  equation    (7)    into   equation    (8)       yields   the 
governing    equation  or   vertical   equation. 


2 

^-^  +    X    ^    N^(z)w^    =    0  (9a) 

,    2  n  " 

dz 


subject  to  boundary  conditions. 


w'  =  0      at  z  =   0  (Sea  Surface) 
w*  =  0      at  z  =  -H  (Sea  Bottom) 


(9b) 
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Equation     (6)    written    in    -che   form, 

7^Pn)t    -    Xnf^^PnH     +     /3(Pn)x     =     «  (10a) 

represen-s    the   Rcssby      wave  equation    foi:    which      solutions   of 
the    form 


p^  =   EXP[i(kx    *-    ly   -   at)]     ,  (lob) 


are  sought  where  k:  and  %.  are  wavenumbers  in  th^  x  and  y 
direction,  respectively  and  a  is  the  frequency.  Substitut- 
ing from  (lOt)  into  (10a)  yields  a  form  of  the  dispersion 
relationship    fcr    baroclinic  Rcssby    waves,    viz. , 


of 


Equation  set  (9)  forms  a  Sturm-Liouville  eigenvalue 
problem   for   the    vertical    motion,      where  X^    are   the    eigen- 

values and  w«  the  correspcnding  eigenf unctions .  A  solution 
is  sought  for  a  finite  number  of  modes,  m.  These  modes  cor- 
respcnd   to    the   first    m  baroclinic   Rossby   wave   modes. 

If  N  (z)  were  independent  of  depth  at  a  fixed  loca- 
tion     and      time,        the      solution    for      w^      could      be      readily 
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determined   analytically      as        tr igonome-cric  functions      of 

(n7rz)/H    (see    for   example    Eedlosky    (1979,    pp.    356-360)),    pro- 
vided 


^n      "   HE    '      ^    "    1>2,3,...     .  (12) 


In  certain  instances  an  analytical  solution  to  equation  set 
(9)  can  be  found  provided  N^z)  is,  cr  approximates,  a  well 
defined  function  (see  for  example  Willmott  and  Mysak(1980) 
who  solve  the  vertical  prcfclera  with  N^z)  as  an  exponential 
function  of  z)  .  Realistic  distributions  of  N'^are  depen- 
dent upon  depth  and  generally  not  well  defined  functions  of 
z,  therefore  the  solution  of  equation  set  (9)  is  obtained 
numerically. 

2.      Determination  of    N"^   and    The    Isotherm    Displacement 

The   Brunt-Vaisala    frequency   squared      is,      as    defined 
earlier,         equal    to  -(g/Po)Pz  •         Since     the    model      neglects 

salinity,    the    simple    state  relation. 


p    =    p^ri   -   a(T   -   T^)]     ,  (13) 


is    used   to    convert  the      Brunt-Vaisala    frequency    squared    from 
a      function    cf      the    vertical      mean    density      gradient    into      a 
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function    of   the    vertical    mean      temperature      gradient,       wher3 
To  is    a      constant    reference  temperatura   and     a.    is   the   coeffi- 
cient of    therial   expansion.  Differentiating   equation    (13) 
with   respect    to    z  and      substituting      inro    the      definition  of 
N^    gives, 

N^    =   ag  T      .  (14) 

^      z 

Equaticn  (14)  is  used  tc  calculate  N^  from  the  mean 
monthly  temperature  fields  output  by  the  ocean  model  for  the 
last  ten  years,  1969-78,  cf  the  model  run.  Only  ^he  temper- 
ature field  in  the  region  18.3  -  56. 9N  and  159.  1E  -  125W 
shown  as  the  enclosed  box  in  Figure  1  is  used  in  the  calcu- 
lation of  N^  AfTier  completion  of  the  numerical  diff  eren-^  ia- 
tion  (right  side  of  equaticn  (14)),  a  ten-year  time  mean,  N-^ 
at  each  grid  pcirt  and  at  19  vertical  levels  (raid-levels  of 
the  20  model  levels)  "  is  calculated.  Finally  N^is  averaged 
over  28  east-west  grid  points,  giving  a  N-^  versus  depth  pro- 
file for  each  latitude  (2C  north-scut.h  grid  points)  .  Justi- 
fication for  the  east-west  averaging  is  based  upon  little 
variability  of  N^ with  longitude  over  the  region  at  a  partic- 
ular   vertical    level.      The    variability   of    N^ with   latitude,    on 
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the  ether  hand,  is  such  that  north-south  averaging  is  not 
justified.  Price  and  Magaard  (1980)  also  found  that  mean  N^ 
profiles  varied  much  more  with  latitude  than  with  longitude, 
and  thus  the  numerical  model  data  are  consistent  with  the 
corresponding  data  taken  from  the  North  Pacific  Ocean.  Pro- 
files of  (N^ r^  ,  averaged  over  longitude,  are  shown  in 
Figures  2  to  5  for  the  latitudes  of  the  grid  points  selected 
in  this  investigation  (triangular  symbols  on  Figure  1).  The 
solid  curves  en  each  figure  are  drawn  to  the  values  of  {t^^^h. 
at    19     vertical    levels      determined    by      aquation    (14).  The 

dashed  curves  are  drawn  to  values  of  (N-^)^''^  calculated  at 
every  five  meters  from  z  =  0  to  z  =  -4000  m  using  interpola- 
tion. The  interpolation  is  done  numerically  by  means  of 
Aitken's  scheme  of  Lagrange  interpolation  described  by 
Hildebrand    (1956,   pp.    49-50)    using   the   values   of    (N''^)^/^at   the 

—  ~   1/2 

19  vertical  levels  and  an  assigned  value  of  (N^ )'  at  -4000  m 
as  tabular  input.  At  latitude  48, 8N  and  54. 8N  th9  assigned 
bottom  value  for  (N^  )^/^  was  given  the  identical  value  calcu- 
lated at  level  19  (2830  m) .  At  the  other  selected  latitudes 
the  assigned  bottom  value  was  determined  by  linear  extrapo- 
lation  using    the    levels    below    1000    m    (levels    15    to    19) . 
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Th€  profiles  shew  a  pronounced  variability  with 
latitude.  Belcw  approximately  900  m,  {^^  y^  is  small  (less 
than  0.001  1/ssc)  and  nearly  constant  (model  simulated 
isothermal  ccndition  in  the  deep  ocean).  At  latitudes  north 
of  U7N  this  constancy  is  evident  at  depths  balow  150  m.  The 
maxiiEum  value  cf  (N-^  )'  (naximum  vertical  -cemperarure  gradi- 
ent) occurs  at  level  2  (16,5  m)  at  latitudes  south  of  22N 
and  between  latitude  40N  and  SON.  At  the  other  la^.itudes 
the  maximuir  occurs  near  the  surface  (5  m)  .  A  secondary  max- 
imum occurs  between  latitudes  30N  and  45N  at  a  depth  between 
200    m   and    400    m    (inferred    location    of    the    mean   thermocline) . 

The  second  field  tc  be  computed  using  the  model  sim- 
ulated mean  monthly  temperature  fields  is  the  isotherm  dis- 
placement. The  isotherm  displacement,  D(2,t)  ,  is  the 
displacement  of  a  monthly  mean  isotherm  from  its  long-term 
average  level.  It  is  defined  mathematically  as  -T'/(T)2 
where  T'  is  the  temperature  fluctuation,  T»  =  T  -  T.  f  and 
(T)^  are  the  long-term  (ten-year)  averaged  temperature  and 
vertical  temperature  gradient,  respectively.  In  the  numeri- 
cal calculation  of  D(z,t)  the  assumption  is  made  that  the 
temperature  field  is  linear  in  z  between  model  levels.  The 
isotherm      displacement   versus      time      curves     at    19      vertical 
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Figure    2.       Picfiles    of    (N^)^at  Laritudes    18. 3N    and   2U.4N 


23 


N    (l/SEC) 

0.005  0.010 


0.015 


0.000 


o 
o 


CL. 
Q 


o 
o 
V  ■ 

(M 

1 


o 
oa 

I 


o 
o 
rvj- 

I 


o 
o 

n 


o 
o 
o 


LiJ 
Q 


LATITUDE:    30.5   N 


o 
o 

U3- 


O 
O 


I 


o 
o 
oo- 

tM 
I 


O 
O 
(M- 
CO 
I 


O 
O 

m 

I 

o 
o 
a 


LRTITUDE:    36.6   N 


Figure  3.   Prcfiles  of  (iJ^  j'^^at  Latitudes  30. 5N  and  36. 6N 
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levels  frcm  tke  ten-year  model  simulation  are  shown  in 
Figures   6,    1,    and  8    for    th€   grid    point   at    30. 5N    139. 1W. 

By  virtue  of  its  definition,  D(z,t)  positive,  in 
general,  inplies  a  temperature  colder  than  the  long-term 
mean.  D(z,t)  negative  implies  a  temperature  warmer  than  the 
long-term      mean.  The    isotherm      displacement      curves,        as 

depicted  in  Figures  6  to  8,  reveal  a  pronounced  seasonal 
signal    (annual    period)         at   levels   above   71.5   m.  At    100   m 

(level  6)  and  telow  the  seasonal  signal  becomes  less  well 
defined  and  losses  definition  with  depth.  The  seasonal  sig- 
nal at  the  upper  levels  is  primarily  the  model  simulated 
mean  monthly  tsmperature  fields  response  to  the  seasonally 
varying  solar  radiation  provided  as  model  input.  Below  100 
m  (Figures  7  and  8)  a  low  frequency  (long  period)  signal  is 
evident.  The    isotherm      displacement   curves      at   the      other 

selected  grid  points  (not  shown)  ,  in  general,  show  similar 
characteristics  of  an  upper  level  annual  cycle  and  low  fre- 
quency signal  at  lower  levels.  The  significant  features, 
proceeding  from  south  to  north  at  longitude  155. 9W,  are 
that,  1)  the  depth  at  which  the  seasonal  signal  remains 
distinguishable  increases  reaching  1760  m  (level  17)  at 
54. 8N,      2)       the    amplitude    of      the   isotherm      displacement      at 
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mid-l-rvels  increases  from  a  maximum  of  approximately  50  m  at 
18. 3N  to    5000    m   at   54. 8N.  This   increase   in   amplitude-    is   a 

result  of  the  near  iscthemal  conditions  evident  below  150  m 
at  the  northern  latitudes  (i.e.  (T)^  is  small  thereby 
increasing   the    value    of    D(2,t)). 

-•      2h^    Solution    of    The   Vertical    Problem 

Equation  set  (9)  is  solved  numerically  using  an 
Interactive  Crdinary  Differential  Equation  (lODE)  solver 
which  uses  a  fcurth-order  Bunge-Kutta  integration  technique. 
The  eigenvalues,  X^ ,  are  found  by  means  of  the  "c-rcot 
method"  descrited  by  Preisendorf er  and  Gonzalez  (1 973)  .  In 
applying  the  "c-rcot  methcd"  equation  (9a)  is  separated  into 
two    first-order    differential   equations. 


Qw^    -   A  (15a) 

dz      ■   ^ 


di      .    _    ^2    ^2f7.w     .•  (15b) 


St      -    -K    N^(z)w^    . 


A   trial      value   is  chocsen    for      \,      and   aquations       (15a)       and 

(15b)      are    integrated  using  lODE    starting    at   the   sea    bottom, 

z    =   -H    (4000    m)  ,      with   boundary    values     of    w^'    =    0   and     0=    1. 

As      the      integration      proceeds  the      zeros      of      the      function 
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function  c  (z)      =    w'(2)       -      0(z)Z(z')       ara    noted.         In    the 

above  expression  Z(z*)  equals  w'(z»)/0(z')  where  z'  is  the 
value  of  z  a-  the  bottom  tcundary  which  changes  as  -he  inte- 
gration continues.  Since  the  boundary  condition  states  that 
w'  (-H)  =  w'(2*)  =  0,  then  2(z)  must  equal  zero  and  the  c 
function  reduces  to  c  (z)  =  w'(z).  The  eigenvalue  problem 
thus    reduces    tc    finding    values  of    z    where    w'=    0. 

For  example,  for  the  (N^^/^  profile  at  30. 5N  and  a 
trial  value  of  X=  2.5  sec/m,  c(z)  =  v*  (z)  =  0  when  z  = 
-2451ffl,  -liCSm,  -607in,  -326m,  and  -112[n.  Varying  X  by  Q.  1 
sec/m,  noting  when  w'=  0,  and  tracing  out  the  points  results 
in  five  curves  (eigencur ves)  as  shown  in  Figure  9.  The 
first  five  eigenvalues  which  satisfy  the  boundary  condi- 
tions, equations  (9b),  are  found  by  the  intersection  of  the 
eigencurves  with  the  z  =  0  line.  The  first  five  eigenvalues 
for  this  example  are  shown  to  the  right  of  the  ordinate  in 
Figure  9;  namely,  X-,  =  0.4735  sec/ip,  X2=  0.8934  s^c/:^, 
X3  =  1.3380  =ec/m,  X^  =  1.7814  sec/m,  X^  =  2.1823  sec/m. 
Justification  for  limiting  the  calculation  cf  eigenvalues  to 
the  first  five  modes  will  be  made  apparent  in  the  next  chap- 
ter. An  advantage  in  finding  the  eigenvalues  in  this  manner 
is    that    the    depth  to    the    sea      bottom   need    not    be    fixed.         If 
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the  d€p-ch  to  the  bottcm  were,  for  example,  3  ka  rather  -han 
U  km,  the  eigenvalues  satisfying  the  same  given  boundary 
conditions  could  be  obtained  by  noting  the  value  of  at  the 
intersection  of  the  eigencurves  and  the  z  =  -1000  m  line. 
The  reciprocal  eigenvalues,  Cp  which  correspond  to  the 
internal  phase  speeds  in  the  absence  of  rotation  associated 
with  each  tarcclinic  mode,  are  given  in  Table  I  for  the 
latitudes  of  each  selected  grid  point.  The  tabular  values 
of  Cn  show  a  decrease  with  latitude  and  wi-h  increasing  mode 
number. 

Osing  expendable  bathythermograph  (XBT)  data  -aken  at 
weather  station  November  (30N,  140W)  during  July  1966  to 
December  1S69,  Emery  and  Magaard  (1976)  obtained  c^  values 
of  c,  =  2.99  ir/sec,  Cj  =  l.'il  m/sec,  c^  =  0.92  m/sec,  C4  = 
0.72  m/sec,  and  Cg  =  0.51  m/sec.  Thess  values  are  slightly 
higher  than  the  values  given  in  Tabl?  I  for  latitude  30. 5N 
(latitude  cf  ncdel  grid  point  closest  to  the  latitude  of 
weather  station  November)  .  The  difference  can  be  attributed 
to  the  deeper  depth  ar  weather  station  November  (4755  m)  and 
to  larger  values  for  N^  (i.e.,  the  values  for  N^calculated 
from  the  model  simulated  temperature  field  is  smaller  than 
observed)  . 


33 
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Figure  9,   Eigenc^orvso  and  Eigsnvaluas  for  Laritude  30. 5N 

Having  determined  the  firsT.  five  eigenvalues  by  -h* 
"c~rcot  method",  the  r.axt  step  is  to  calculaxe  the  ccr- 
respcading    eigenf unctions,     w^   ,    using   lODE.     The 
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TABLE  I 
Reclprccal  Eigenvalues  (in/=9c)  for  Modes  1  to  5 


LATITUDE 

1 

MODE 
2      " 

NUMBER 

T" 

4 

5 

18.3N 

2.2578 

1.2442 

0.8505 

0.6274 

0. 50  4  9 

2U.4N 

2 .  23  3  1 

1.1966 

0.8202 

0.6099 

0.4922 

30.  5N 

2.11  19 

1.  1193 

0.7474 

0.5614 

0.4582 

36. 6N 

1.7937 

C.9774 

0.6245 

0. 4719 

0.3883 

4  2.7N 

1 .  32  8  9 

C.7510 

0.4736 

0.3620 

0. 3077 

48.  8» 

0.9299 

0.48  31 

0.3294 

0.2629 

0.2264 

54.  8N 

0  .  55  1 1 

0.2598 

0.2153 

0. 1722 

0.  1356 

eigenfunct icns   are  then   ncrmalizGd,    following   Emary  and 
Magaard  (1976),  as 


Y„(2) 


w^'(z) 


(16) 


f      N^(z)(w^'(z))^dz]^ 


/2 


-H 

where  n  represents  tha  mode  number.  Since  equation  S'rt  (9) 
represents  a  Sturm-Licuville  system,  the  eigenfunct  ions  are 
orthogonal  with  respect  to  the  weight  function  N^  (z)  over 
the    vertical   cclumn    z   =    -H   to   z    =0,    viz. , 


/ 


N^(z)w    •(z)w    '(z)dz    =    0    ;         n    i'   m 
n  m 


(17) 


-H 
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I 

Tharefore,      ncrmalization      of    the   eigenf anc-ions      rssal-s    in 


the  crthonormal  sst^Jy^l  ,  viz. , 


/ 


N      n     m  nm 


-H 


where  6nm  ^^  ^^®  Kronscker  delta  function.  A  more  detailed 
discussion  ccnceining  the  orthogonaiity  condition  for  tha 
eigenfunctions  is  given  by  Willmott  and  Mysak  (1980) .  The 
normalized  eigenfunctions  for  the  seiectad  latitudes  are 
shown  in  Figures  10  through  16.  The  characteristics  of  the 
eigenfunctions  are  related  to  the  profile  of  N^  (z) .  At 
latitudes  south  of  UON  (Figures  10  to  13)  the  nodal  points 
of  the  normalized  eigenfunctions  are  concentrated  near  the 
upper  half  of  the  column  in  response  to  the  model  simulated 
upper  ocean  N^ maximum  (maximum  temperature  gradient).  At 
the  northern  latitudes  (Figures  15  and  16),  the  normalized 
eigenfunctions  are  very  nearly  sinusoidal,  particularly  at 
the  higher  [redes.  This  corresponds  to  the  near  constant 
profile  of  (N^  j^^  (Figures  4  and  5)  below  150  m. 

The  analysis  procedure  is  concluded  by  fitting  the 
normalized  eigenfunctions  to  the  simulated  isotherm  dis- 
placement,  D(z,t).     This  is   done  in   the  time   domain  by 
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NORMRLIZED  VERTICAL  VEL0CITY(5EC/Mxxl/2 ] 
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Mode  1:  solid;  Mode  2:  Dashed;  Mcde^S:  Long  dash/Short  dash; 
Mode  H:    DoTt€d;   Mode  5:  Long  dash/Dot 

Figure  10.   Ncrmalized  Sigenf unc tions  at  18. 3N 
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NORMALIZED  VERTICAL  VEL0CITY[SEC/Mxxl/2 ) 

-20.0     -10.0       0.0       10.0       20.0 


Mode  1:  Solid;  Mode  2:  Dashed;  Mode  3:  Long  dash/Short  dash; 
Mode  U:  Dotted;   Mode  5:  Long  dash/Dot 

Figure  11.   Ncrmalized  Eigenf unctions  at  24. 4N 
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NORMRLIZED  VERTICRL  VEL0CITY(SEC/Mxxl/2) 

-20.0     -10.0       0.0       10.0       20.0 


Mode  1:  solid;  Mode  2:  Dashed;  Mode  3:  Long  dash/Short  dash; 
Mode  4;  Dotted;   Mode  5:  Long  dash/Dot 

Figure  12.   Normalized  Eigenf unctions  a^  30. 5N 
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NORMALIZED  VERTICRL  VEL0CITY(SEC/Mxxl/2 ) 

-30.0    -20. 0       -10.0  0.0  10.0         20.0         30.0 


Mods    1:    Solid:    Mode    2:    Dashed;    Mode    3:    Long    dash/Short    dash; 
Mode  4:     Dotted;      Mode    5:    Long    dash/Dot 

Figure    13.      Ncrmalized  Eigenf unc tions   at    36.  6N 
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NORMALIZED  VERTICAL  VELOCITY  (  5EC/Mx><l/2  ] 
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noae     :3uxiu^^^^  ^^  Dotted;   Mode  5:  Long  dash/Dot 

Figure  14.   Normalized  Eigenf unctions  at  n2.1^ 
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NORfinLIZED   VERTICRL   VELOCITY  [5EC/Mxxl/2) 
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Mode    1:    solid;    Wcde    2:    Dashed;    Mode    3:    Long    Jash/Shor-    dash; 
Mode  4:     Dotted;      Mode    5:    Long    dash/Dot 

Figure    15.      Ncrmalized  Eigenf unc tions   at    48. 8N 
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NORMnLIZED  VERTICRL  VE:L0CITY[5EC/Mxx1/2  J 
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Mode  1:  Sclid;  Mcde  2:  Dashed;  Mods  3:  Long  dash/Short  dash; 
Mode  4:  Dotted;   Mode  5:  Long  dash/Dot 

Figure  16.   Normalized  Eigenf unctions  at  5a. 8N 
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making  us€  cf  the  ccmple-eness  of  the  eigenf unctions  10 
9quat6  D(z,t)  to  the  sum  of  the  normalized  eigenf  unctions 
multiplied    by    an    amplitude   coefficient   A^,    Z^z.  , 


D(z,t)    =   y^       A^(t)Yj^(z)     . 


(19) 


n=l 


The    coefficients,   A^  (t) ,      are   readily   obtained   by    making    use 
of    the   ort hcncrmali ty  condition.  Multiplying    each    side    of 

equation     (19)     ty    N^ymSna    then   integrating    from   z   =   -H    t.o    z   = 
0    results    in 


"       ■  09 

I     N^(2)Yjjj(z)D(z,t)dz    =  V"        [A,^(t)      r    N^(z)Yj^(z)Yjj^(3)dzl 
-H  "=^  -H 


(20) 


Application  cf  the  othogorality   condition  (18)   ro  equatio; 
(20)  yields 


/ 


A  (t)  =   /   N^(z)y  (z)D(z,t)dz  .        (21) 
n        I         n 


-H 


Equaticn  (21)  is  \:s€d  to  ccmpute  A^(t)  for  thr  first 
five  modes  for  each  selected  grid  poin^.  The  devia"*-icn  from 
the    icdai    fit,    as  defined    fcy    Emery    and   Kagaard    (19  76)  ,  eqn.(22), 
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is  also  calculated.  The  deviaticr.  -cerm  qudnzi^axiv^sly  shows 
how  well  the  model  simulated  isotherm  displacement  flucrua- 
tions  are  fit  by  the  first  five  baroclinic  modes.  A  value 
of  DEV(t)  =  0  implies  that  the  iso-cherm  displacement  fluctu- 
ations can  he  fully  explained  by  the  first  five  baroclinic 
modes. 


o 


/' 


2(z)[D(z,t)    -   V^'  A    (t)Y    (z)]^    dz 
/  ^      n  n 


„                                           n=l 
DEVCt)    =    — .  (22) 


/ 


-H 


N^(z)    D^(z,t)dz 
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III.       RESULTS 

The  coefficients.  An  (t) ,  and  the  devia-ion,  DEV(t)  ,  cal- 
culated at  each  selected  grid  point  are  shown  in  Figures  17 
to  2^.  At  lower  latitudes  (Figures  17  to  19  and  Figure  24) 
rhe  aitplitude  of  coefficient  A,  is  generally  larger  than  the 
coefficients  A2  tc  A5  indicating  the  importance  of  -he  first 
mode   at    the   southern    latitudes.  At   latitude    24, UN    (Figure 

18)  a  seasonal  signal  is  evident  for  coefficients  A^  and  A5 . 
The  seasonal  signal  of  the  higher  modes  becomes  more  appar- 
ent with  increasing  latitude  indicating  that  the  seasonally 
varying  model  input  parameters  appear  to  be  reflected  in  the 
higher  modes  at  higher  latitudes.  At  the  northern  most  grid 
point,  54. eN  (Figure  23)  the  overall  amplitude  increases  an 
order  of  magnitude.  Additionally,  the  pattern  of  each  coef- 
ficient appears  nearly  identical  with  only  slight  differ- 
ences  in    magnitude. 

The  deviation  term  at  latitudes  18. 3N  and  24. 4N  remains 
generally  belcw  50%  during  most  of  the  ten-year  model  run, 
although  sciie  tiire  periods  are  higher.  North  of  24. 4N  there 
is      a     general   increase      over      the     ten-year   period      in      the 
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Figure    17.      A^  (t)    and   DEV  (t)     at    18.  3N    155.  9W 
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Figure    18.       k^{t)    and   DEV  (t)    at    24. UN    155. 9W 
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Figure    19.       A^  (t)    and   DEV  (t)     ar    30. 5N    155.  9W 
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Figure  21.   A^,  <t)  and  DEV(t)  at  42. 7N  155.  9W 
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Figure  22.       A_(t)    and   DEV  (t)     at    U8.8N    155.  9W 
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Figure   24.      A^  (t)    and  DEV(t)    a-    30. 5N    139.  IW 
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I 

deviation    curv<=    with    the    exception      of    the    two    northern    ^cst 

grid  points  (Figures  22  and  23)  .  As  a  comparison,  the  first 
five  coefficients  and  the  deviation  term  computed  by  3mery 
and  Kagaard  (1976)  using  U8  months  of  data  for  weather  sta- 
tion Novemter  are  shown  in  Figure  25.  The  results  shown  in 
Figures  24  (model  grid  point  closest  to  weather  station 
November)  and  25  compare  reasonably  well.  Both  show  the 
largest   amplitude      of  coefficeint      \  .  The    U8-raonth      mean 

deviation  computed  for  weather  station  November  was  32% 
while  the  ten-year  mean  deviation  is  38.2^  for  model  grid 
point    30. 5N    139. 1W. 

To   examine  the  relative  contribution      made   by    each    baro- 
clinic   mode      to    fittirg    the   simulated      isotherm   displacement 
fluctuations       (equation  (1 9) )  ,      the     ten-year    mean      deviation 
(equation     (22))       was    calculated    for    each    selected    grid    point 
adding  one    mode   at  a    time.  The  results    are   given  in   Table 

II.  Column      1       (under    modes)  gives      the      ten- year      mean 

deviation  when  only  the  first  baroclinic  mode  is  fit  to  the 
isotherm      displacement    fluctuations.  Cclamn      2    gives      the 

ten-year  mean  deviation  when  modes  1  and  2  are  fitted  to  the 
isotherm  displacement  fluctuations  and  so  on.  The  tabular 
values  show      the   importance  of  mode      1   at    the      selected   grid 
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Figure  25.      An  (t)  ,   and  DEVift)     for   weather    sta 
(from  Emery    and   Magaard,    1976) 


tion   November 


points  south  cf  30. 5N  (i.e.,  over  40%  of  the  deviation  is 
accounted  for  ty  the  first  baroclinic  moda  alone) .  The  ccn- 
tribution  mad*  by  the  fifth  barcclinic  mode  in  reducing  the 
mean  deviation  is  less  rhan  6%  at  the  three  southern  ip.cst 
grid  points,  and,  except  for  the  northern  most  grid  point, 
is  less  than  IC^  elsewhere.  Based  upon  these  results,  it  is 
apparent  that  adding  higher  baroclinic  modes  (modes  greater 
than  5)  wculd  add  little  tc  reducing  the  deviation  term  com- 
pared to    mcde   number    one. 
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TABLE    II 
Mean    Eeviaticn     (?)    at   Selected  Grid   Poin-cs 


GRID    POINT 


18. 3N    155. 9W 


24. UN    155. 9W 


30. 5N    155. 9W 


36. 6N    155. 9W 


4  2.7N    155-9W 


4  8. 8N    15  5.9W 


54. EN    155. 9W 


30.  5N    139. 1W 


MODES 
1  1+2         U2+3 

56.3  46.0  36.4 

54.1  45.9  39.0 


U2-^3  +  4 
28.6 
35.0 


1+2+3+4+5 
26.8 
31.6 


57.4 

54.2 

50.9 

46.  9 

41.0 

65.2 

55.7 

53.3 

45.3 

38.8 

96.2 

68.1 

60.7 

54.2 

46.2 

56.  1 

34.5 

28.8 

24.  1 

21.6 

90.  1 

57.8 

43.0 

26.6 

10.0 

63.6 


i6.5  50.9 


41.  1 


38.2 


The  mean  deviation  (column  5  of  Table  II)  generally 
increases  with  latitude  up  to  42. 7N  but  then  decreases  at 
the  two  highest  latitudes.  Additionally,  the  contribution 
of  modes  2  to  5  in  reducing  the  mean  deviation  is  signifi- 
cant at  the  ncithern  most  points.  To  examine  the  relevance 
of  these  results,  the  cutoff  period  (minimum  period)  for  the 
propagation  cf  free  Rosshy  waves  of  mode  n  (see  Emery  and 
aagaard  (1976))  given  by  T  =  (47rfAn)//3  was  evaluated.  The 
results  are  shewn  in  Figure  26.  The  dotted  line  delineates 
the    annual    period  and  the    dashed    line   the    five    year   period. 
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Figure   26.      Cutoff  Period    versus    La-itude 
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The    curves    cf   Figure    26    shew   the    following    results: 

•  Eased  upon  the  theoretical  cutoff  oeriod  anriual  neriod 
tarocli.iic  free  Rossby  waves  may  sxis-c  nor-h  cf  25n  but 
only    mode    1    can    occur. 

•  Earoclinic  free  Rossby  waves  at  -che  annual  period  can- 
not 


were    prmarily    evicl 

in    the    North  pacific   Ccean)  . 

•      Earoclinic    free      Rossby   waves    at      the   five      year    period 
cannct    exist  north    cf  approximately   52N. 

The  results  given  above  are  significant.  Since  the  free 
Rossty  wave  analysis  was  conducted  using  only  ten  years  of 
model  simulated  temperature  data,  free  Rossby  waves  of  peri- 
ods greater  than  approximately  five  years  can  not  be  ade- 
quately resolved.  Therefore,  the  results  obtained  for  model 
grid  point  54. EN  155. 9W  losses  significance.  Likewise,  at 
grid  point  48. 8N  155. 9W  the  theoretical  cutoff  period  curves 
show  that  only  first  mode  taroclinic  Rossby  waves  at  periods 
greater  -^han  approximately  three  years  can  exist.  There- 
fore, the  low  deviation  at  latitude  43. 3N  resulting  prima- 
rily from  modes  2  to  5  is  not  necessarily  due  to  baroclinic 
Rossty  waves.  The  opposite  may  be  said  of  the  selected  grid 
points  at  latitudes  13. 3N  to  36. 6N.  At  these  latitudes  the 
first  baroclinic  mode  is  dominant  and  the  cutoff  period  is 
4.5  months  at  18. 3N  and  12.8  months  at  latitude  36. 6N.  At 
latitude    42. 7N   the   cutoff    period      is    21.4    months.         However, 
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only  the  first  and  second  modes  can  exist  at  periods  less 
than  five  years.  The  mean  deviation  is  high  when  only  the 
first  and  second  modes  are  firtsd  to  rhe  isotherm  displace- 
ment fluctuations  (96.2%  and  88.15S  in  columns  1  and  2, 
respectively,  of  Table  II).  Thus  at  model  grid  point  42. 7N 
155. 9W  the  model  simulated  isotherm  displacemen-c  fluctua- 
tions cannot  be  explained  by  free  baroclinic  Rossby  waves  in 
the   time    dcmain. 

Thus  far  the  examination  of  baroclinic  free  Rossby  wave 
signature  in  the  model  siiulatad  temperature  field  (or  iso- 
therm displacement)  has  teen  limited  -co  the  time  domain 
only.  In  order  to  conducx  an  analysis  of  baroclinic  free 
Rossty  wave  signature  in  the  space  domain  requires  carrying 
out  a  cross-sprsctral  fit  cf  the  x  and  y  direction  wave  num- 
bers at  the  appropriate  frequencies.  In  the  context  of  this 
study  a  cursory  examination  of  the  wave  number  vectors  will 
be  made.  The  locus  of  wave  number  vectors  for  mode  1  at 
appropriate  frequencies  for  latitudes  18. 3N  and  36. 6N  are 
shown  in  Figures  27  and  28.  The  locus  of  wave  numbers  was 
determined  using  equation  (11)  for  three  different  frequen- 
cies. 
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Figure   27.      Lccus  of    Have    Number    Vectors    for    Mode    1    at    18. 3N 
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Figure   28.      Lccus  of    Wave    Number    Vectors    for    .lode    1   at    36. 6N 
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The  spacing  between  model  grid  points  u-cilized  in  this 
study  is  cr  the  order  of  226  km  in  the  north/south  dirsc-icn 
and  251  km  (at  36N)  in  the  east/west  direction.  The  shor- 
test wavelength  resolvable  by  the  model  is  equal  to  twice 
the  spacing  between  grid  poin-s ,  which  in  this  case  is  on 
the  crder  of  approximately  500  km.  Therefore,  Rossby  waves 
resolvable  by  the  model  must  have  wavelengths  grea-^er  than 
or  equal  to  the  2^6.  value,  where  Ad  is  the  grid  spacing. 
The  wave  numbers  in  the  x  and  y  direction  which  correspond 
to  a  wavelength  which  is  greater  than  tne  2Ad  wavelength  are 
located  withir  the  dotted  rectangle  in  Figures  27  and  23. 
Stated  ancther  way,  wave  numbers  outside  of  the  boxed-in 
region  correspond  to  wavelengths  too  small  to  be  resolved  by 
the   mcdei. 

The  inner  circle  in  Figure  27  is  the  locus  of  wave  num- 
bers found  whan  a  frequency  corresponding  to  a  six-month 
period  is  us  =  d.  The  middle  circle  and  ou-cer  circle  are  for 
frequencies  corresponding  to  periods  of  nine  months  and  12 
months,  respectively.  Note  that  due  to  -he  different  scale 
of  eaca  axis  the  circles  appear  as  ellipses.  Portions  of 
all  three  circles  fall  within  the  enclosed  region  indicating 
that    wavelengths   greater      than   the    2Ad      wavelength   can    occur 
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at  th€  af  or6m€r.t  Icned  frequencies  (periods).  The  circles  in 
Figure  28  are  icr  frequencies  corresponding  to  periods  of  18 
months  (inner  circle),  24  months  (middle  circle)  and  30 
months  (cuter  circle) ,  respectively.  It  is  noted  that  at 
36. 6N  only  waves  having  a  period  of  about  30  months  or 
greater  have  wavelengths  long  enough  to  be  resolvable  by  the 
model.  It  fcllcws  that  waves  of  44d  or  84d  would  have  even 
greater  periods.  Therefore,  increasing  latitude  requires  an 
increase  in  the  Rcssby  wave  period  to  allow  for  wavelengths 
long  enough  tc  be  resolved  by  the  model  grid  spacing.  Addi- 
tionally, using  equation  (11),  it  may  be  shown  that  for 
wavelengths  that  are  resolved  by  the  model  grid,  k^  +  e^«  Xp  f^ 
and  thus  the  period  2  tt /a  must  be  larger  than  the  cutcff 
period.  At  high  latitudes  such  long  periods  would  not  be 
well  resolvad  ty  the  ten-year  simulation.  Thus,  even  though 
low  deviations  result  at  high  latitudes  (column  5  of  Table 
II)  the  isotherm  displacement  variability  cannot  be  attrib- 
uted   to    actual   Rossby   waves   with    any    degree   of   certainty. 
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IV.       CONCLU SIGNS    AND    RECOMMBN D ATIQN S 

Eased    upcn   the  preceeding   analysis   the    following   results 
are    summarized: 

•  On  the  average,  free  baroclinic  Rossby  wave  signature 
in  th€  jcdel  simulated  isotherm  disolacement  fluctua- 
tions is  evident  in  the  time  domain* for  the  selected 
grid    pcin-s    south  of    UON. 

•  Az  the  selected  grid  points  south  of  40N  the  first 
haroclinic    node    shows   predominance. 

•  Eased  upon  the  locus  of  wave  number  ellipses,  Rossby 
waves  at  the  appropriate  periods  are  long  enough  to  be 
resolvable    witnii:  the   model    grid   spacing. 

•  Eased  upon  the  deviation  term^  certain  times  within  the 
ten-V€ai:  nodel  run  period  indicate  that  the  first  five 
baroclinic  modes  do  not  adequately  fit  the  isotherm 
displacement  fluctuations  (i.e.,  the  deviation  is  high, 
therefore,  the  isotherm  displacement  fluctuations  can- 
not te  fullv  explained  by  baroclinic  free  Rossby 
waves) . 

•  At  selected  grid  pcints  north  of  40N  the  analysis 
results    are    inconclusive. 

•  There  is  reasonably  gccd  agreement  between  the  analysis 
results  at  model  grid  point  30. 5N  139. 1W  and  that  con- 
ducted by  Emery  and  Magaard  (1976)  for  weather  station 
November     (3CN    140H)  . 

The    results   stated   in    point   number    five   above,    are   a    product 

of    several   factors.         The   model    run   time    (ten    years)       sets   a 

constraint    upcn   the    periods   that      may    be    adequately   resolved 

by    the   analysis    method.         The   cutoff   periods    of    the   selected 

grid      points    north      of      ION  gradually      begin      to   approach      a 

five-year    period,   and   wavelengths   resolved   by   the   model   grid 
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result  in  periods  considerably  greater  than  the  cu-off  (min- 
imum) period.  Periods  greater  than  five  years  can  nc-^.  be 
adequately  resolved  by  the  analys:.s  method.  A  longer  run 
time  is  required.  Secondly,  since  the  ocean  model  is  a 
"fresh"  water  model  the  effect  is  to  dimenish  the  magnitude 
of  the  Brunt  —  Vaisala  frequency,  pariiicularly  '.z  higher 
latitudes.  This  results  in  eigenvalues  that  are  large  com- 
pared to  "real"  ocean  values  (  \  is  inversely  proporricnal  to 
N)  .  The  larger  eigenvalues  then,  in  turn,  increase  the 
value  of  the  cutoff  pericd  more  than  whaz  occur  when  only 
the  Coriolis  parameter  is  considered.  Finally,  a  seasonal 
signal,  resulting  from  the  seasonally  varying  model  input 
parameters,  dcminates  the  isotherm  displacement  fluctuations 
over  a  greater  portion  of  the  water  column  at  the  northern 
selected  grid  points.  Thus,  the  seasonal  signal,  which  is 
primarily  a  thermodynamic  effect,  contaminates  the  first 
five    modes. 

Eased      upon   these      results      several  recommendations      for 
further   research    may    be    set  forth: 

•  The  seasonal  signal  resulting  from  thermodynamic 
effects  should  be  removed.  This  could  be  achieved  by 
either  neglecting  the  model  layers  above  the  main  ther- 
mccline  or  by  applying  a  filter  technique  to  remov=  the 
the    regular     seasonal  signal.  The    former      suggestion 

would  be  impractical  at  the  high  latitude  model  grid 
points,  however,  due  to  the  penetration  of  the  seasonal 
signal    down    to    the    bottom  of   the    water   column. 
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•  This  study  set  the  m-ran  curren-  equal  to  zero.  Further 
studi€S  should  include  the  mem  current  in  +he  analy- 
sis, particularly  in  the  region  of  the  Califcrnia  cur- 
rent, and  North  Pacific  Current.  A  study  by  Kang  and 
Maaaard  (1979)  has  shown  that  the  North  Pacific  Current 
has  a  significant  imfact  on  the  first  baroclinic  shear 
mode . 

•  Since  the  model  neglects  salinity  further  research 
should  study  the  feasability  cf  incorporating  salinity. 
This  would  provide  for  more  realistic  density  profiles 
and  improve  the  analysis  results  at  t ne  higher 
latitudes . 

•  The  mechanisms  which  produce  Eossby  waves  in  the  North 
Pacific  Ocean  is  still  not  fully  understood.  The  anal- 
ysis has  shown  that  certain  time  periods  within  the 
trodal  run  period  south  of  latitude  i+ON  cannot  be 
exTslained  fully  by  the  first  five  baroclinic  modes. 
Apparently  some  other  factor  is  involved.  Further 
research  could  be  done  to  isolate  any  peculiarities  in 
the  the  synoptic  wind  field  used  as  model  input  since 
the  six-hourly  synoptic  wind  was  the  only*  varying 
parameter  from  one  year  to  the  next.  The  study  should 
also  te  extended  to  investigate  the  source  or  Rossby 
waves  (i.e.,  whether  they  result  from  large-amplitude 
variations  of  the  wind  stress  curl  near  -tae  eastern 
boundary  as  hypothesized  by  White  and  Saur  (1981)  or 
whether   seme  oth«r    mechanism   is   involved). 


•      All    further    studies    should    include  comparison   to    appro- 
€      observational     studies        wh 
by   include    a   model   evaluation. 


priate      observational     studies        when    "available,    ""and 
thereby   i]     '    " 


This  study  has  been  a  first  step  in  analyzing  the  low  fre- 
quency dynamics  of  the  Haney(1980)  ocean  model  below  the 
mixed  layer.  A  complete  understanding  of  the  dynamics  of 
the  large-scale  variability  of  the  ocean  temperature  field 
below  the  oixed  layer  (thermocl ine)  is  the  primary  goal. 
This  understanding  will  in  turn  lead  to  improvements  in  the 
ocean  model  and  its  diagnostic  and  prognostic  capabilities, 
the    ultimate   gcal. 
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